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Abstract 

We report on a calculation of e' k /ek at next-to-leading in the 1/N C expan- 
sion and to lowest order in Chiral Perturbation Theory. We discuss the 
short-distance matching, the scale and scheme dependence as well as the 
long-distance short-distance matching. We include the two known chiral 
corrections to our result and discuss further order p 4 corrections to it. 
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1 Introduction 



Recently, direct CP violation in the Kaon system has been unambiguously estab- 
lished by the KTeV |]T] experiment at Fermilab and by the NA48 experiment 
at CERN. Their results together with the previous NA31 || and E731 ||] mea- 
surements produce the present world average 



Re {s'k/sk) 



;i9.3±2.4) • 10' 



(1) 



Further reduction of the statistical error to the order of 1 ■ 10~ 4 is expected. 

A lot of effort has been put in the theoretical side to get a Standard Model 
prediction for this quantity in the last twenty five years. Recent reviews and 
predictions are || |6|, 0, [| |9], [10], [Ll|. Here, we would like to report on a calculation 
of this quantity in the chiral limit and next-to-leading (NLO) order in 1/N C [Q. 
We also discuss how it changes when the known chiral corrections, i.e. final state 
interactions (FSI) and tt — 1] mixing are included. Comments on the different 
approaches to obtain the contributions from Q G and Q 8 to e' K will also be given 
in the Summary. 

Direct CP-violation in K — > titt decays amplitudes is parameterized by 
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A [K L -> (7T7T) J=2 ] 
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A [K S - (7T7T) 7=0 ] 

We want here to predict K — > at NLO in 1/N C and to lowest order in the 
chiral expansion. 

In the isospin symmetry limit, K — > tttt invariant amplitudes can be decom- 
posed into definite isospin quantum numbers as [A = —iT] 
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with 5q and 62 the FSI phases. 

To lowest order in Chiral Perturbation Theory (CHPT), i.e. order e°p 2 and 
e 2 p°, strong and electromagnetic interactions between n, K, 77 and external 
sources are described by 



£( 2 ) = -JLtr (uX + X +) + e 2 C 2 tr (QUQU+) 



(4) 



with U 



2 = exp(A a 7r a /F ) and = iv) \D^U)v) . \ a are the Gell-Mann matri- 
ces and the ir a are the pseudoscalar-mesons tt, K, and r\. Q — diag(2/3, —1/3, —1/3) 
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is the light-quark-charge matrix and x+ = u^x u ^ J r u X^ u an d X = 2-Bo diag(m„, ma, m s ) 
collects the light-quark masses. To this order, F n = F = 87 MeV is the pion 
decay coupling constant. Introductions to CHPT can be found in Refs. fl"3|, [TJj. 



To the same order in CHPT, the chiral Lagrangian describing \AS\ = 1 is 



*-|AS|=l 



+CF* 



CF$e 2 G E tr (a 32 Q) 
G 8 tr (A 32 uX) + G^tr (A 32 ^ 



+G 27 f j * l tr (AijU M ) tr (A kl u^ 



+ h.c. 



(5) 



with A i:j = u\iju\ (Xij) a b = S ia 8 jb , Q = u ] Q 



C 



3G, 



v ud v: s « 



-1.06 • 1(T 6 GeV~ 



5^2 

The SU(3) x SU(3) tensor t 1 ^ can be found in [15]. Using this Lagrangian, 



(6) 
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(7) 



and 80 = 62 = 0. In the presence of CP-violation, the couplings G$, G27, and G E 
get an imaginary part. In the Standard Model, ImG 2 7 vanishes and ImGg and 
lmG E are proportional to Imr with r = —\ t /\ u and Aj = Vi^V* s and where V^- 
are Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. 



2 Short-Distance Scheme and Scale Dependence 

Strangeness changing transitions happen in the Standard Model by the exchange 
of one W^-boson. This fact implies that all physics between and 00 has to be 
taken into account when calculating Kaon to pions amplitudes. In particular, 
this implies the intervention of strong interactions at all energies, i.e. from the 
perturbative region to the non-perturbative one. 

The two very different scales involved in K — > tttt, i.e. the iy-mass and the 
Kaon mass make effective field theory methods very useful. This is standard and 
we outline the steps needed. Below the IV-mass, the effective action Tas^i is 
obtained by integrating out the heavy particles, top, Z, and iy-bosons and using 
the operator product expansion (OPE). The leading contributions consist of four- 
quark operators and higher dimensional operators are suppressed by M^. Using 
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the renormalization group equations, Pas=i is brought down to some perturbative 
scale below the charm quark mass. Perturbative matching and OPE is used at 
thresholds of the successive heavy particles which are integrated out. This full 
process implies several choices of short- distance schemes, regulators, and operator 
basis. Of course, physical matrix elements cannot depend on these choices. 

Explicit calculations have been performed including gluonic and electromag- 
netic Penguin-diagrams to one-loop in [|l^] and to two-loops in |l7j in two schemes. 



These are the naive dimensional regularization (NDR) scheme and the 't Hooft- 
Veltman (HV) scheme. The choice of the operator basis and other choices like 
infrared regulators can be found in those references. 

The Standard Model Ta5=i effective action at scales v below the charm quark 
mass, is 



18 
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Ta5=i ~ fd*xQ i (x)+h.c. (8) 

i=i 

where Cj = Z{ + Tjji are Wilson coefficients and Qi{x) are four-quark operators. 
At this level, we have resummed the large logarithms [a s {p) \og{v / Mw)} n and 
a s [ce a (v) log(v/Mw)) n to all orders in n. 

At low energies, it is more convenient to describe the AS = 1 transitions with 
an effective action T^ =1 which uses hadrons, constituent quarks, or other objects 
to describe the relevant degrees of freedom. A regularization scheme like an 
Euclidean cut-off which separates long-distance physics from the short-distance 
physics which is integrated out and working in four dimensions is also more 
practical, as well as another operator basis, for instance the color singlet Fierzed 
one. The effective action Ta5=i depends on all these choices and in particular 
on the scale \i c introduced to regulate the divergences generated analogously as 
(H) depends on the scale v. It also depends on effective couplings analogously 
to the Wilson coefficients in (|8|). As usual, matching conditions have to be set 
between the two effective field theories. This is done by requiring that 5*-matrix 
elements of asymptotic states are the same at some perturbative scale 

(2|r£2 =1 |i) = (2|r AS=1 |i>. (9) 

Both sides are separately scale and scheme independent and therefore the short- 
distance scale and scheme dependences are consistently treated. This matching 
is done at the perturbative level using the OPE in QCD. Notice that even if the 
regulator chosen is the same in both sides, there can be finite terms appearing 
in the matching. The matching conditions fix analytically the short-distance 
behavior of the couplings g>j 

&Gu c ,---) = .f(c^), «»,-•■)• (io) 



This was done explicitly in |T9| for AS* = 2 transitions and used in [12| for 
AS" = 1 transitions. To the best of our knowledge, these are the first places were 
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the short- distance scheme dependence is analytically cancelled in the final matrix 
element at next-to- leading order in the 1/N C expansion. Notice that in the whole 
process we never required the matching of scales in two different regularization 
schemes as is sometimes stated in the literature for this type of 1/N C calculations. 

2.1 The heavy X-bosons Method 

We find it convenient to use an effective field theory of heavy color-singlet X- 



bosons coupled to QCD currents and densitities as argued in [19], |20], [21], f22| . 
One advantage of this is that two-quark currents are unambiguously identified 
and that QCD densities are much easier to match than four-quark operators. E. 
g., the operator 

Qi(x) = [s7 M (l - 7 5 )d] (x) [u7^(l - 7 8 )w] 0) 
is reproduced by the effective action 

Tx = ^(/i c ,---)/d 4 l/Xf{[^(l-75)^(x) 

+ [u 7m (1- 75 )u]0e)} (11) 

We use an Euclidean cut-off at \i c in four dimensions as regulator of the UV 
divergences, with this the X-boson effective action is completely specificied. This 
regulator is very convenient to separate long- from short-distance physics. In (|TTD , 
the higher than \x c degrees of freedom of quarks and gluons have been integrated 
out. 

The complete X-boson effective action for the Standard Model AS = 1 tran- 



sitions can be found in [12]. We are now ready to calculate consistently AS 



Green functions with the X-boson effective theory. 



3 Long-Distance— Short-Distance Matching 

We study the two-point AS = 1 Green function 

U{q 2 ) =i J d 4 xe iqx (0\T{Pl{0)P J {x)e irx \0) 

where are pseudoscalar sources with the quantum numbers to describe K — ► 7r 
amplitudes. 

After reducing the external legs one gets K —>■ it off-shell amplitudes. Taylor 
expanding them in external momentum and 7r, and K masses one can obtain the 



couplings of the CHPT Lagrangian [15, 23]. Then one predicts K — > nn at a 



given order. This procedure is unambiguous. One can use also K — > n and K — > 
vacuum transitions form lattice simulations to predict the lowest order CHPT 
couplings |TR EI El . 
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Up to know we have not used any 1/N C argument. We count the gi couplings 
as effectively order 1 since they contain the large logarithms of the running be- 
tween My/ and fi c even though some of them start at NLO order in 1/N C . At 
leading order in the 1/N C expansion the contribution to the Green's function (pj) 
is factorizable. It only involves two-point functions. At this order the result is 
model independent and the scale and scheme dependence is exactly taken into 
account. However, we know that the factorizable contribution fails to reproduce 
Rea by a factor around 6 and Rea 2 by a factor around 1/2. 

At the same order, Q 6 contributes to G% proportional to 

|^| 2 ~C 6 (^)(0|gg|0) 2 (z/)L 5 (/i) (12) 
and Q 8 to e 2 G E proportional to 

|2 8 | 2 ~C 8 (z/)(0|gg|0)V). (13) 
To leading order in 1/N C , the scale dependence of the Wilson coefficients cancels 



exactly the one from the quark condensate pq , p6| . The scale fi dependence of 
L5 requires NLO in 1/N C contributions to cancel. In fact, at this order there 
also appears an IR divergence in the factorizable contribution of Qq to G$ that 



necessarily requires non-factorizable contributions to cancel |22[ . 

The non-factorizable contribution is NLO in the 1/N C expansion and involves 
the integration of four-point functions Up.p.j a j b over the X-boson momentum 
r E that flows through the currents/densities J a and J&. It can be schematically 
written as 

n(g 2 )~ / ^n PiPjJaJb (q E ,r E ) . (14) 

We separate long- from short-distance physics with an Euclidean cut-off fi in r E . 
The short distance part of the integral from /i up to 00 can be treated at NLO 
within OPE QCD using a cut-off regulator if fi is high enough. 

It has been emphasized that dimension eight operators may be numerically 
important when the cut-off scale fi is of the order of 1 GeV El], |27j . This issue can 



be treated in a straightforward way in our approach. We have to use the OPE in 
the short-distance part up to the required dimensions. In fact, those contributions 
are under control at NLO in 1/N C using factorized matrix elements. We stopped 
at dimension six operators for the present results, we plan to check the effect of 
dimension eight operators elsewhere. 

Up to here, there is no model dependence in our evaluation of K — > it ampli- 
tudes at NLO in 1/N C within QCD. 

What remains is the long distance part from up to fi. For very small values 
of fi one can use CHPT and still the result is model independent. However, 
CHPT at order p 4 starts not to be enough already at relatively small values 
of fi, i.e. ~ 400 MeV. To match with the short- distance part, we need clearly 
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to go beyond this energy. The first step is to use a good hadronic model for 



intermediate energies. The model we used is the ENJL model described in [28 



It has several good features -it includes CHPT to order p , for instance- and also 
some drawbacks as explained in JH|. Work is in progress to implement the large 



N c constraints along the lines of |29j . There is very interesting work on a different 
approach to the calculation of weak matrix elements within the 1/N C expansion 
in [|2l], [3C|| . It has started to produce the first results for three- and four-point 
Green's functions. In particular the determination of the parameter in the 
chiral limit is clean |31[ . Their result agrees well with the chiral limit calculation 
in 



19, 21] done within our present approach. 



in the Chiral Limit 



To a very good approximation [14, 18 



■K\ 



1 Rea 2 

y^Reao I Rea ' Rea 2 



Ima Ima2 



(15) 



The isospin amplitudes ai are defined in ([3]). The lowest order CHPT values for 
Re a and Re a 2 can be obtained from a fit to K — > mr and K — > nuir amplitudes 
up to order p 4 in [^ . 

Our results in [22| reproduce the value for Reao finding a nice matching of 
short-and long-distances. For the coupling G27, which dominates Re 02, we don't 
find such good stability though the behaviour is much better than the quadratic 
divergence found in || [33j . In Figure [l] we give G 2 7 and in Figure [| we present 
ReG 8 . 

The AI = 1/2 enhancement is reproduced within 40 %. However, due to 
the present lack of a complete understanding of Rea 2 , we prefer to use the ex- 



perimental values for the lowest order CHPT amplitudes Re a/ [32) to predict 



■k- 



At lowest order in CHPT, the imaginary part of Gg is almost all from the 
Qe operator with very small contributions from Q3, Q4, and Q5. Then to a very 
good approximation 



ImG« 



80 t n 
-— lm.Ty 6 (v) 
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L 5 (fl) B 6 (fi, v). 



(16) 



At large iV c , B^fL^v) = 1. At the same order, the imaginary part of c 2 Ge is 
almost all from the Qg operator with a very small contribution from Qj. Then 
to a very good approximation 



Im (e 2 G E ) 



-5Imr?/ 8 (z/) 
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Figure 1: We plot the result for G27. Curves I and II are different ways of doing 
the two-loop a s running and mul and add are the resummation on the Wilson 
coefficients is done either multiplicatively or additively. The curve quadratic is 
the result of references || [£| . 



fe|0») 



B*(u) 



(17) 



and B 8 (u) = 1 at large N c . With the present values of the CKM matrix elements 



Imr 



-6.72 ■ 10" 



We use the chiral limit MS quark-condensate |3j, 

(0|gg|0)(lGeV) 



£0(1 GeV) = 



Fo 2 



;i.75±0.40) GeV. 



(19) 



Its leading scale dependence is analytically canceled by the one in the Wilson 
coefficients Cq and C% pB, 26[. The uncertainty induced in e' K by B is around 40 



% and often not included in error estimates. Incidentally, higher order corrections 
to the contribution of Qq to G$ do not change the chiral limit condensate into 
the strange quark condensate as is commonly used ||22|| . 

There is no NLO in 1/N C correction to c 2 Ge from Q 8 in nonet symmetry. 
This is a model independent result. In the real world, the U(1)a anomaly gives 
a mass to the pseudoscalar singlet field and octet symmetry is a good symmetry. 
Therefore, there are 1/N% corrections to this result which we intend to study 
within our l/N c approach elsewhere BH] using data to obtain B$ pi EU1 BBI as 
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Figure 2: We plot the result for ReGs- The labels are like in Figure [p. 



well as a lowest meson dominance model. They will give the size of the U(1)a 

studied 



anomaly corrections to our NLO in 1/N C results. The results in [ |10l , 
the B 8 parameter with the additional assumption that dimension six operators 
dominate the OPE. 

Our values of the relevant bag parameters in the NDR scheme at 2 GeV are 
in Table |l|. 



Table 1: Results for B 6 {u)^ NDR , B 7 {p)W NDR and B 8 (u)^ NDR at v 
2 GeV. The WS value m s (2GeV) 
rescaling the results in [pBl |37|. 



(119 ± 12) MeV was used in fTTJ and for 



(1/2)NDR 
6 



(3/2)NDR 
7 



(3/2)NDR 



Method 



B 



(2GeV) B 



(2GeV) B. 



(2GeV) 



NLO 1/N e [x limit] g 
LMD |3 



Dispersive [^6 



QCD Sum Rules 
CHPT NLO 1/N C [| 
Lattice |38fl 
Chiral Quark Model 
FSI Omnes |1] 



2.5 ±0.4 



1.0 ±0.5 
1.5 ~ 1.7 



1.2 - 
1.55 



1.7 



0.8 ±0.1 
0.94 

0.78 ±0.17 
0.7 ±0.2 
-0.1 ~ 0.09 
0.5 ~ 0.8 
~ 0.9 



1.35 ±0.20 

1.6 ±0.4 
1.70 ±0.39 
0.4 ~ 0.7 
0.7 ~ 1.1 
~ 0.9 
0.93 



Details and the input parameters of the determination of e' k /ek are in [ |I~2|j . 
Here we only give the results and main conclusions. To lowest order in 1/N C and 
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in the chiral limit, we get 



' K 



Ek 



(33 - 16) • 1(T 4 = (17 ± 7 ± 6) • 10" 
(17 ±9) • 10~ 4 



(20) 



The first error is from the uncertainty in the value of the quark condensate and 
the second is a typical 1/3 on the factorizable contributions. We cannot give here 
an error estimate due to the non-inclusion of non-factorizable 1/N C corrections 
since they are a new type of contributions. 

Including our calculated NLO in 1/N C non-factorizable contributions we get 



= (83 - 23) • 10~ 4 = (60 ± 24 ± 20) • 10~ 4 
= (60 ± 30) • 10" 4 



(21) 



again in the chiral limit. The first error is once more due to the uncertainty in 
the quark condensate. The second one is an estimate of the model uncertainty 
i.e. around 1/3 of the factorizable contribution plus 40 % of the non-factorizable 
one added quadratically. Notice that there is a factor larger than three between 
( 20|) and (21]). This result is plotted in Figure where one can appreciate the 



quality of the matching. 
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Figure 3: We plot our result for ImGs at NLO in 1/N C . The labels of the curves 
are as for Figure 0. 



To get the result in (|2T| ) we have used our chiral limit NLO in 1/N C determi- 
nations of Im Gs which we show in Figure |] and of Im G e which is in Figure [5] 
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Figure 4: We plot our result for ImGs at NLO in 1/N C . The labels of the curves 
are as for Figure |2|. 



5 Higher Order CHPT Corrections 

The role of final state interactions in the standard |5], [J predictions of s' k /sk has 
been recently clarified in |39|] . See also f4Cj . 

We have taken a different strategy to predict e' K . We start with Equation 
QlBI) . The ratio Ima//Rea/ has no final state interactions to all orders. One 
of the main problems at present of K — > 7r7r lattice calculations is that they 
cannot include FSI. There are very interesting recent developments towards a 



full calculation within lattice QCD including FSI, see and references therein 
for previous work. 

It would be very interesting to see how much of the ratios 1m ai/Ke a j could 
still be obtained from lattice simulations without FSI. One could use CHPT to 
eliminate analytically the FSI contributions for both the real and imaginary parts. 

The only place where FSI are present is in the normalization factor Re /Re do, 
but one can take as a first step its experimentally known value. In fact, lattice 
calculations can provide us with an estimate of s' k /ek to lowest order CHPT 
with K — > 7r and K — > vacuum results fl2f . This would be very interesting. Of 
course, one also wants to predict the AI = 1/2 enhancement ratio Re 02 /Re ao 
but K — > 7T calculations look more straightforward in the lattice at present ||24| . 

About isospin breaking due to m u 7^ m^, only effects due to tt° — 77 mixing 



2^, |42] are under control. They have recently calculated to CHPT order p |l2] . 
Though tt°—i] mixing is all from m u 7^ isospin breaking at lowest CHPT order, 
this is not true at higher orders. In fact, due to the smallness of the order p 4 term 
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Figure 5: We plot our result for ImGg at NLO in 1/N C . The labels of the curves 
are as for Figure |2|. 



in the n° — rj mixing contribution, it could be that there are other m u — rrid effects 
equally important at the same order, so that a full p A calculation is mandatory 
This effect adds a contribution to Im a 2 which is parameterized usually 



as 



[Im d2 



IB 



n 



Re Gt2 

0.16 ±0.03 m 



IB 



Im do 
Re do 



(22) 



We will use the value CIjb 

Other purely real p 4 corrections even in the isospin limit are mostly unknown. 
They have been taken into account within the approach in M and partially in 



Electromagnetic corrections have been considered in p6| |44| . Very little is 
known at present of the size of the order p 4 contributions and beyond. Obviously 
there is more work to be done in this direction. 

To lowest order in 1/N C but including FSI and tt° — rj chiral corrections -notice 
that these two corrections are actually NLO in 1/N C - we get 



(24-15.5) • 10~ 4 
(8.5 ±4.6) ■ 10~ 4 . 



(8.5 ±3.5 ±3) • 10" 



(23) 



The errors correspond to the same discussion of the chiral limit results ([ZD) 
and in particular do not include the error of the non-included non-factorizable 
contributions. This result is model independent. 
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Notice that small changes in the isospin zero or in the isospin two contribu- 
tions are extremely amplified, so that a 20 % correction to both but in opposite 
directions can translate into a factor two in the value of e' K . 

Including our calculated NLO in 1/N C non-factorizable contributions, we get 



= (59.6 - 25.6) • 1(T 4 = (34 ± 14 ± 11) • KT 4 
= (34 ±18) ■ 10~ 4 . 



(24) 



Notice that the two known chiral corrections reduce the value of e\ 



K • 



6 Summary 

We have reported on a calculation of s' k /ek at next-to- leading in a l/N c ex- 
pansion and in the chiral limit ||12|| . Emphasis has been given to discuss the 
short- distance scheme and scale dependence matching in Section |2| and to the 
long-distance short-distance matching in Section [| 

We have also given the two model independent results at leading order in the 
1/N C expansion and in the chiral limit in (|20|) and including the two known chiral 
corrections, namely FSI and 7r° — 77 mixing in (23). 



We also have given the main result of our work which is the value of e' k /sk in 
the chiral limit and at next-to- leading order in a 1/N C expansion (|21~|) . Including 
again the two known higher order CHPT corrections to this result, we obtain 
our final number for e' k /ek in (|24]) . Notice however, as discussed in Section 
|^, that most of the higher order CHPT corrections are actually unknown and 
as remarked in several recent work they could affect numerically the Standard 
Model value of e' k /e k significantly fl2], |43], |44 |. 



We would like to discuss the determinations of the two main bag parameters 
for e' k in Table [l]. About B & , which is the most uncertain, at present, there are 
only two estimates of the non-factorizable corrections not proportional to quark 
masses to this important parameter. They are the results in and in [ |12|| . 

The Dortmund group II has included non-factorizable corrections to B e in 
the chiral limit. For that, they use CHPT to order p 4 and find a quadratic 
dependence in the cut-off scale. They have also included partially higher order 
CHPT corrections to that result. 



To the best of our knowledge, the result in ||12|1 , is the first one that includes 
NLO in 1/N C non-factorizable corrections to B 6 in the chiral limit and find a 
logarithmic matching with the Wilson coefficients see Figure |3] for e' k /ek and 
Figure |] for the imaginary part of G%. 

The Trieste group || uses B e = 1 in the chiral limit and all the difference 
between the factorizable result and their final number quoted in the Table [I] 
are higher order CHPT corrections. They have not calculated non-factorizable 
corrections in the chiral limit. 
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Notice that the calculation of Bq in [|10j only includes the factorizable contri- 
butions. 



In [11], the leading order in 1/N C values of Bq = 1 and B% — 1 are corrected 
with 7T7T final state interactions using the Omnes solution. In our opinion, this 
is at present the best way of taking into account within CHPT those important 
corrections. 

For Q 8 we go at NLO in 1/N C but in the chiral limit, and we find that, as 



explained in Section |] and ]T2[ . 

Bg DR (2GeV) x = 1.35 ± 0.20 (25) 

is a model independent result. There are however 1/N% corrections induced by 
the U(1)a anomaly. These effects were estimated in Jl0|, |36| with the additional 
assumption that dimension six operators dominate the OPE in this case. We 
plan to investigate this issue within our approach [[35| . 

Again, only NLO corrections proportional to quark masses are included in 



the result of the Trieste group to B 8 ||. The results from the lattice |38| do not 
include FSI but we don't known how much do they include of other higher order 
chiral corrections. 

Let us conclude that a lot of advances have been done towards the prediction 
of e' k /ek within the Standard Model. For us, there are two main questions: 
which is the value of B 6 in the chiral limit and which are the non-FSI chiral 
corrections to it. This work tries to answer the first of these questions. There is 
going on a lot of activity to clarify these issues by several groups and we foresee 
a very exciting nearby future. 
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